Light-emitting current of electrically driven single- phot on sources 



O . David M.-T. Kuo 

O 

Department of Electrical Engineering, National Central University, Chung-Li, Taiwan, 320, Republic of 
. ! China 



-a 

G 

O 

o 



X 



(February 6, 2008) 



Abstract 



Time-dependent tunnelling current arising from the electron-hole recombination of 
the exciton state is theoretically studied using the nonequilibrium Green's function 
technique and the Anderson model with two energy levels. The charge conservation 
and gauge invariance are satisfied in the tunnelling current. Apart from the classical 
capacitive charging and discharging behavior, interesting oscillations superimpose on 
| the tunnelling current for the applied rectangular pulse voltages. 
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I. INTRODUCTION 



Recently, much effort has been devoted to the study of single-photon sources (SPSs) made from 
quantum dots (QDs), which have potential applications in quantum cryptography and quantum 
computing 1-3 ). The second order correlation function of light emitted from such devices must be 
measured to determine their functionality as SPSs, that is, the light source antibunching feature must 
be examined. The antibunching features of SPSs were demonstrated in Ref. [1,2], where electrons and 
holes in the QD are excited by optical pumping. Nevertheless, only a few experiments have employed 
electrical pumping is to demonstrate the antibunching behavior of SPSs. From a practical point of 
view, it is more efficient to use the electrically driven SPSs. A prototype SPS made from individual 
QDs embedded in a semiconductor p-n junction was suggested by Imamoglu and Yamamoto 4 . In 
addition to antibunching, enhancing the spontaneous emission rate is necessary, because it is crucial 
to SPS performance. 

For small size QDs, the strong three dimensional confinement effect creates large energy separations 
among the low-lying confined levels. Consequently, it is possible to inject electrons and holes into 
the corresponding ground state energy levels of a single QD and generate photons via electron-hole 
recombination in the exciton, positive /negative trion or biexciton states in the QD. To fabricate 
a single of nanometer-sized QD at a specific location is one of the most challenging processes in 
the realization of electrically driven SPSs. Self-assembled quantum dots (SAQDs) combined with a 
selective formation method may overcome such a difficulty 5 . Even though, some theoretical studies 
have been devoted to electrically driven SPS devices, 6 " 8 the dynamic properties of an electrically 
driven SPS are still not clear. In ref. 4, the Monte Carlo method was used to simulate junction 
dynamics, where time-dependent junction voltage is calculated. Other studies emphasized steady 
state characteristics such as the exciton assistant tunnelling current 6 , Purcell effect on the tunnelling 
current 7 and electrode effects on the exciton complexes 8 . 

The main purpose of this study is to investigate the time- dependent tunnelling current (or sponta- 
neous emission rate) arising from electron-hole recombination of exciton state by using nonequilibrium 
Green's function technique 9-11 , which has been applied to several different systems. 12-15 . This study 
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attempts to clarify how the carrier tunnelling process, applied modulation voltage and temperature 
influence the spontaneous emission rate in an electrically driven SPS. The schematic band diagram 
describing the studied system is shown in Fig. 1. Only the ground states of conduction and valence 
band of the QD are considered. The energy level of E 2 is 20 meV below the Fermi energy level of 
left electrode E F>e = 50 meV. On the other hand, the exciton state denoted by E ex is 15 meV above 
the Fermi energy of right electrode Ep^ = 50 meV. r 2 (i) denotes the tunnelling rate from the QD to 
the left (right) electrode. The emitted photons will be observed when the periodic square bumplike 
modulation v(t) added into the right electrode injects holes into the exciton resonant level. 

II. HAMILTONIAN 

An Anderson model with two energy levels in terms of electron picture is used to describe the 
system as shown in Fig. 1, 

H = e k,La{a k + £k,Rblh + e i4 d i (!) 

k k i=l,2 

+ \i 2 d\d 2 + X^d^di + 5^ifc,i&fcdi + h.c 

k,l 

+ E'l,2°l rf 2 + h - c 

k,2 

where a\{ak) and b\{bk) create (destroy) an electron of momentum k in the left and right system 
electrodes, respectively. The free electron model is used in the electrodes in which electrons have 
frequency-dependent energies e^L = — cj/2 and €k,R = Sk + oj /2 + v{t). Time-dependent modulation 
v(t) denotes the time-dependent applied voltage in the right electrode. d\ (di) creates (destroys) an 
electron inside the QD with orbital energy q = E^ — (— l) l u/2. In this study % — 1 and % = 2 represent, 
respectively, the ground states of valence band and conduction band of individual QDs. The fourth 
and fifth terms describe the coupling of the QD with electromagnetic field of frequency uj. A = — fi r S 
is the Rabi frequency, where /j, r =< f\r\i > is the matrix element for the optical transition and S 
is the electric field per photon. t fc i describes the coupling between the band states of electrodes and 
energy levels of the QD. Note that a unitary transformation, S(t) = exp^nCLM^-^+dl^-d^) ? 
has been used to obtain Eq. (1) via 
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H = S^HmS -iS^—S. 

v 1 dt 

For small semiconductor QDs, the particle correlation is strong. This indicates that the intralevel 
Coulomb interactions of the QD, Uu and U 22 , can not be ignored. In order to simplify problem, it is 
assumed that a regime of the applied voltage can not overcome the charging energies U n and U 22 . s 
To investigate the exciton assistant process 6 ' 8 , the interlevel Coulomb interaction Ui 2 is taken into 
account in Eq. (1) 

H v = U 12 d\dJ 2 d 2 , (2) 
which is invariant under unitary transformation. 



III. FORMALISM 

Although a spontaneous emission process of photons is the quantum effect of an electromagnetic 
field 8 , the electromagnetic field is still treated as a semiclassical field in Eq. (1). The approach 
detailed in ref.fll] is employed to study the spontaneous emission process of photons. The optical 
susceptibility of QDs not only provides the absorption coefficient and the refractive index, but also 
determines the spontaneous emission rate. 11 Therefore, the optical susceptibility of individual QDs 
is calculated using the lesser Green's function defined as Gfj(t,t') = i < d\(t')dj(t) >. Based on 
Keldysh's equation, we have 

G<(M) =/dti/ cft 2 GT iB (t,t 1 )E< m (t 1> t 2 )G^,.(t 2 ,0, (3) 

where £< m , G\ n and G^ ■ are the lesser self-energy, the retarded Green's function and the advanced 
Green's function, respectively. The Einstein summation index is used in Eq. (3). The spontaneous 
emission rate will be suppressed if tunnelling rates T 2 and r\ are smaller than spontaneous emission 
rate R eh . The detailed expression of R e h is given later. To avoid the suppression of the spontaneous 
emission rate, the device shown in Fig. 1 favors the condition of T 2 > F 1 » R eh or F 1 > T 2 » R e h 
to function as a light emitting source. The condition of T 2 > Ti >> R eh , \j will be regarded as a 
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small parameter in the comparison with t k>i . Consequently, the lesser self-energy of Eq. (3) is mainly 
attributed to tunnelling 

Ef i (t,0 = </^r i / i (c) e - fc («-«'), (4) 

where the tunnelling rates are Tj = J2ktk,itt ti S(e — and /j(e) = \/{e e ~ IJ ' F ' i — 1) is the Fermi 
distribution function of the electrodes, in which the chemical potential is given by jiF,i = Ef,% — (l) l u/2. 
It is worth noting that energy-independent tunnelling rates are assumed in Eq. (4) for the sake of 
simplicity. Inserting Eq. (4) into Eq. (3), we obtain two off-diagonal lesser Green's functions 

/df 
^{rifi(e)A r 21 (e, t)A a n (e,t) (5) 

+ T 2 f 2 (e)A r 22 (e,t)A a 21 (e,t)}, 

and 

GUt,t) = i J ^-{r 2 f 2 (e)A{ 2 (e,t)AUe,t) (6) 
+ r 1 f 1 (e)A r 11 (e,t)A a 12 (e,t)}, 

as well as two diagonal lesser Green's functions 

Gf A (t,t) = i J ^{rx/x^fc ^(M) (7) 
+ r 2 f 2 (e)A r 12 (e,t)A a 21 (e,t)}, 

and 

-{T 2 f 2 (e)A r 22 (e,t)A a 22 (e,t) (8) 
+ r 1 / 1 (e)^ 1 (e,*K 2 (e>*)}> 

where J dhe^'^G^h) = A^t) and /dtie-^-^G^Mi) = A ^(e,t). The expression of 
Alj a \e, t) depends on the detailed form of retarded and advanced Green's functions. Eqs. (5) and 
(6) determine the optical susceptibility of QDs. Eqs. (7) and (8) determine the electron occupation 
numbers. Once G 21 (t,t) is determined, Gf 2 (t,t) is also obtained as a result of G 21 (t,t) = (Gf 2 (t,t)Y. 
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To solve G 21 (t,t), some approximations are considered in the following derivation due to the com- 
plicated Hamiltonian of Eq. (1). Because the energy level of i=2 emerges into the band width of 
the left electrode and T 2 is larger than r 1; the retarded and advanced Green's functions for i=2 can 
be regarded as a steady state solution. That is G r 22 a \t,t') = q= i9(t - t') e - i ( e iTr 2 /2)(t-t') ) w hi c h gives 
A r 22 a \e, t) = l/(e - e 2 ± iT 2 /2). As for G r ff \t, t'), the off-diagonal Green's function (i ^ j) are solved 
using Dyson's equation, 

Gf\t,t>) = J dt 1 G%\t,t 1 )\ iJ G%\t 1 ,t'). (9) 

Due to Xij < tk,i and Glj a \t, t') in terms of the first order parameter A, the Green's functions G^ 
in Eq. (9) are A-independent functions. Consequently, A r 2 ^\e, t) = X 2 \A 22 l \e)A r 1 [ a \e,t) and Eq. (5) 
is rewritten as 

Gii(M) = a 2 i / ^{mrMU^A^e) (10) 

+ f 2 (e)T 2 \Al 2 (e)\ 2 A a u (e,t)}. 

To obtain the spontaneous rate resulting from the imaginary part of G 21 (t,t), the imaginary part on 
the both sides of Eq. (10) is used 

ImG^t) 

= A21 / ^{/iWril^Ccor^^W -A a 22 (e)} (11) 

+ / 2 (e)r 2 |^ 2 ( e )| 2 i[^ 1 (e,t)-A r n (e,t)]}. 
Employing G\ — G\ — Gf — Gf, where the greater Green's function is given by 

G?j(M) = JdhJ dt 2 GlSt^)K, m {tiMG a m ,{t 2 ,t) (12) 
with the greater self-energy 

^(t,0=^/^r i (l-/ l ( e ))e-^'), (13) 

Eq. (11) can be simplified as 



ImG< h (t,t) (14) 

= -Aefc / ^{(i -/ e (c) -A( e ))r,|^( e ,t)| 2 r e |^ e ( e )| 2 }. 

In Eq. (14) the electron and hole picture are employed to label i=2 and i=l, respectively. Note 
that Eq. (14) contains a factor of (1 — f e (e) — A(e)). It is always possible to write \ e hImG< h (t, t) = 
X e (uj, t) — X a (uj, t), where X a (u, t) and X e (u, t) are, respectively, in proportion to (1 — / e (e))(l — /h(e)) 
and / e (e)/h(e). X a (u,t) and X e (u,t) denote the absorption and emission spectra, respectively. We 
only focus on the emission spectrum, which was frequently reported in experiments, 

Xefat) = K h j ^{/e(6)/,(6)r,|^(6,t)| 2 r e |A: e (6)| 2 }. 

By reference to ref. 11, the current arising from the electron-hole recombination of the exciton state 
may be written as 

J sp (t) (15) 
= eaJdco co 3 J ^mf h (e)T h \Al h (e,t)\ 2 Im[A a e M 

with a = 4n;!/i 2 /(6c 3 ft 3 eo), where n r and e are the refractive index and static dielectric constant of 
system, respectively, e denotes the electron charge. Note that Eq. (15) uses r e |Ag e (e)| 2 = 2 Im[A® e (e)] 
which is valid only for the steady state. It should be noted that A r hh (e,t) satisfies the condition of 
T h \A r hh (e, t)\ 2 = —2 Im[A r hh (e,t)\ — d ^ Ahh ^^ resulting from total charge conservation and gauge 
invariance 16 . p(uj) = oj 2 of Eq. (15) arises from the density of states of photons. J sp (t) is determined 
by the time-dependent interband joint density of states and the factors of f e (e)fh(e). According to 
the expressions of X e (u, t) and J sp , we prove that the intensity of emission spectrum is linear variation 
with respect to J sp . 3 

To solve the spectral function of holes A r hh (e,t), the retarded Green's function of holes G r hh (t,ti) 
is derived to obtain 

G r hh (t, h) = (1 - N e )g r h (e h , t, h) + N e g r h {t h + U eh , t, h) (16) 

with 
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g h (eh,t,ti) = -i6(t-ti)e Jt i , (17) 

where eh = —E h + ^, and N e is the electron occupation number at steady state. Two branches exist in 
Eq. (16) ; one corresponds to the resonant energy level of E h with a weight of (1 — N e ), and the other 
corresponds to the exciton resonant level of E ex = E h — U eh with a weight of N e . Consequently, holes 
injected into the energy levels of QDs depend not only on the applied voltage of the right electrode, 
but also on the electron occupation number N e , which is given by N e = J ^■f e (e)Im e _ E ^_ iT j 2 . As 
long as the applied voltage is insufficient to inject holes into the resonant energy level of Eh, it is only 
necessary to consider the exciton branch of N e g r h (th + U e h,t,ti) in Eq. (16). The detailed expression 
of G hh (t,ti) depends on the applied voltage v(t), considering a rectangular pulse with the duration 
time of As and amplitude A as v(t). For t < t < t + As, according to Eq. (17), 

Afct) (18) 
N e A 



1 + 



e-e h + iF h /2 K e-e h - A + iT h /2 

(1 _ e i(.e-e h -A+ir h /2)(t-t )^y 



For t > As, 



A h (e,t) (19) 

:i + 



e-e h -A + iT h /2 x t-e h + iT h l2 

(I _ e i(e-e h +ir h /2)(t-As)^ _ e i(e-e h ~A+iF h /2)Asyy 

The time- dependent tunnelling current associated with the spontaneous radiative transition in indi- 
vidual quantum dots is next discussed based on Eqs. (15), (18) and (19). 

IV. RESULTS AND DISCUSSION 

Due to current conservation, the time-dependent tunnelling current is given by J(t) = J sp (t) = 
eT s (t), where T s (t) denotes the time-dependent spontaneous emission rate. Owing to small tunnelling 
rates (Y e /EF, e « 1 and T^/Ep^ « 1), Eq. (15) can be reduced to 



J sp (t) = 2e*R eh f e (T)A h (T,t), 



(20) 



with 



Mt)= I -fh(e)r h \A r h (e,t)\ 2 

J 71 



(21) 



In Eq. (20) we define the time-independent spontaneous emission rate R e h = a Ql x , where Q ex = 
Eg + E e + E h — U e h- According to Eq. (20), the time-dependent feature of J sp (t) is determined by 
the hole occupation number of Eq. (21). 

To reveal the time- dependent behavior of J sp (t), Eq. (20) is solved numerically, and shows J sp (t) 
for various duration times at zero temperature in Fig. 2 for the tunnelling rate Th = 0.5 meV and 
A = 20 meV; the solid line (As = to), the dashed line (As = 3 to) and the dotted line (As = 5 to). 
The current almost reaches saturation at t — 5 to- Apart from the classical capacitive charging and 
discharging behavior (exponential growth and decay) 17 , the interesting oscillations superimpose on 
the tunnelling current not only for t < As but also for t > As. However, the amplitude of oscillation 
for t < As is small. In particular, the oscillation period for t > As does not depend on the magnitude 
of duration time. The oscillatory current is yielded from the particle coherent tunnelling between the 
electrodes and the QD. Such a coherent tunnelling was also pointed out as the mechanism of quantum 
ringing for electrons tunnelling current through the single dot embedded in an n-i-n structure 18 . In 
ref. 18, oscillatory current is not observed as t > As. According to the result of quantum interference 
between the outgoing wave (leaving the QDs) and the wave reflected from the barrier, this oscillatory 
behavior could be observed in the discharging process 19 . 

To examine if the oscillation currents shown in Fig. 2 depend on applied bias strength, Figure 
3 illustrates J sp (t) for different applied biases with duration time As = 3 to at zero temperature 
and T h = 0.5 meV. The solid line denotes A = 15 meV, and the dashed and dotted lines denote, 
A = 20 meV and A = 25 meV, respectively. For amplitude A = 15 meV, the Fermi energy of the 
right electrode just reaches the alignment with the resonant exciton level. J sp (t) displays a strong 
oscillation characteristic. If the pulse amplitude is increased to A = 20meV, the resonant exciton 
level is covered by the right electrode reservoir, and therefore the magnitude of J sp (t) increases. When 
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the pulse amplitude to A = 25 meV, J sp (t) becomes saturated . The results shown in Fig. 3 clearly 
indicate that the oscillation feature of J sp (t) depends on the magnitude of A. Due to the emitted 
photon numbers in proportion to J sp (t), the features of emitted photon numbers as functions of time 
will correspond to the results shown in Fig. 3. 

In Figs. (2) and (3) hole tunnelling rate is set as T h = 0.5 meV, it is interesting to understand 
how the tunnelling rate influences J sp {t). J sp (t) for different hole tunnelling rates is shown in Fig. 4; 
the solid line (T h = 0.5 meV), the dashed line (T h = 0.75 meV) and the dotted \me(T h = 1 meV). 
Increasing the tunnelling rate, J sp {t) reaches exponential growth quickly. This also indicates that 
electrically driven SPSs can quickly reach the maximum photon emission efficiency within a shorter 
operating time compared to optically driven SPSs with a phonon bottleneck. To clarify, the relation 
between the photon number N s and J sp {t) should be constructed. Because dN ^ = r sp (t), the time- 
dependent photon number N s (t) = J^dt T sp (t) is obtained. For the step-like pulse (As = oo,), 
A r h (e, t) is given by Eq. (18). Due to the small oscillatory amplitude in the charging process, Eq. (18) 
is approximated as A h (t) = f h (T)(l - exp~ Tht ) for t > 0. Consequently, N s (t) = R eh f e (T)f h (T)t(l + 
(exp- Fht - l)/(iy)). We found that N s (t)/t = R eh f e (T)f h (T) as T h t » 1. This indicates that a 
higher tunnelling rate allows an electrically driven SPS to quickly reach the maximum photon emission 
efficiency. 

Finally, the finite temperature effect on J sp (t) is shown in Fig. 5; the solid line (ZcbT = 0), the 
dashed line (ksT = 1 meV), the dotted line (ksT = 2 meV) and the dot-dashed line (ksT = 3 meV). 
Fig. 5 shows that the spontaneous photon emission rate T sp (t) is suppressed by temperature effects 
resulting from a factor of f e (T)f h (T) and iV e 2 (T) (see Eqs. (18) and (19)). iV e 2 (T) is yielded by the 
electron-hole interaction. From Fig. 2 to Fig. 5, J sp (t) is in units of J = 2e R e h- For InAs QDs, the 
typical value of R eh is on the order of 1 [ieV (I / ns) . 20 This indicates that J sp results from spontaneous 
photon emission on the order of nA, which can be readily measured. 
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V. SUMMARY 



The expression of tunnelling current J sp (t) arising from spontaneous radiative transition in in- 
dividual quantum dots is obtained by using the nonequilibrium Green's function method. J sp (t) is 
found to be functions of the spontaneous emission rate R e h, the tunnelling rates T h (T e ) and the fac- 
tor of f e (T)f h (T). In addition to the exponential growth and decay features corresponding to the 
charging and discharging process, the oscillatory behavior of J sp {t) as a function of time is observed 
as a result of particle coherent tunnelling between the electrodes and QDs, which is suppressed by 
the temperature effects. 

In this study, resonant tunnelling carriers are used to yield the triggered single-photon sources, 
in contrast to the captured carriers typically used in an optically driven SPS. 21 ' 22 Due to the phonon 
bottleneck effect 23 , it is predicted that the QD capture rate of electrons will be low. This could 
reduce the performance of SPS devices, which use captured carriers as the source of single-photon 
generations. Using carriers injected via the resonant tunnelling process can prevent such problems. 
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Figure Captions 

Fig. 1. The schematic band diagram for the single quantum dot embedded in a p-i-n junction. 
The ground state energy level of electrons is 20 meV below the Fermi-energy of left-electrode, and 
the exciton state for holes is 15 meV above the Fermi-energy of right-electrode. The periodic square 
bumplike modulation v(t) supplies holes into the quantum dot. 

Fig. 2. Time-dependent current arising from the electron-hole recombination J sp (t) for different 
duration times of applied rectangular pulse voltage with amplitude A = 20 meV at zero temperature. 
Current is given in units of J = 2e R e h, where R e h denotes the time and temperature-independent 
spontaneous emission rate, and time is in units of to = h/meV. 

Fig. 3. Time-dependent current arising from the electron-hole recombination J sp (t) for different 
applied voltages with duration time As = 3 t at zero temperature. Current is given in units of 
J = 2e R e h, where R e h denotes the time and temperature- independent spontaneous emission rate, 
and time is in units of t = h/meV. 

Fig. 4. Time-dependent current arising from the electron-hole recombination J sp (t) for different 
tunnelling rates of holes at zero temperature and the applied voltage with duration time As = 3 t 
and amplitude A = 20 meV. Current is given in units of J = 2e R e h, where R e h denotes the time 
and temperature-independent spontaneous emission rate, and time is in units of to = h/meV. 

Fig. 5. Time-dependent current arising from the electron-hole recombination J sp (t) for different 
temperatures at the applied voltage with duration time As = 3 to and amplitude A = 20 meV. 
Current is given in units of J = 2e R e h, where R e h denotes the time and temperature-independent 
spontaneous emission rate, and time is in units of t = H/meV. 
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